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Abstract 



(<— ^ i Results of some recent investigations of democratic three-body decays of 

t^^ I nuclear systems are discussed. In particular, we consider experimental studies 

^-p ■ of three-body {a + N + A^)-decay of j4 = 6 nuclei and J-matrix calculations 

Qv , of monopole excitations in ^^C and El transitions in ^^Li in the three-body- 

continuum cluster models a + a + a and ^Li -\- n + n, respectively. 

^' 1. Introduction 



Investigation of the decay properties of nuclear resonant states is an important 
part of the nuclear structure studies. Usually in such investigations only two-body 



Vh ' decay channels are allowed for. Three-body decays are much more informative than 

the two-body ones, but the full description of a three-body decay is a very compli- 
cated problem. Thus it is very important to derive simple but accurate approximate 
approaches for the analysis of the experimental data and for making reliable theo- 
retical predictions. One of the approaches of such kind is the approach based on the 
democratic decay approximation. 

A three-body continuum spectrum wave function in the asymptotic region is gen- 
erally a superposition of the components describing two- and three-body decay chan- 
nels ^. The democratic decay approximation accounts for the three-body channel 
only; two-body channels associated with the appearance of bound two-body sub- 
systems and 'non- democratic' subdivision of the system of the type (2+1), are not 
allowed for in the approximation. Therefore, this approximation is valid only for 
the study of a three-body system in the case when all two-body decay channels are 
closed and the only open channel is a three-body channel. Generally speaking, if 
the three-body channel is the only open one, than the system can have two types of 
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asymptotics ^. One of the asymptotics corresponds to the situation when one of the 
particles is scattered by another one and the third particle is a spectator. From the 
general physical point of view, this asymptotics is supposed to be of little importance 
for a nuclear system excited in some nuclear reaction and decaying via a three-body 
channel. The alternative type of the asymptotics is a superposition of ingoing and 
outgoing six-dimensional spherical waves (or, equivalently, six-dimensional plane wave 
and ingoing or outgoing six- dimensional spherical wave) and corresponds to the situ- 
ation when the decaying system emits (or/and absorbs) three particles from one and 
the same point in space. The version of the scattering theory that allows only for 
three-body asymptotics of this type is called true three-body scattering theory, or, 
alternatively, 3 — *> 3 scattering theory; the approximation based on the allowance for 
3 — !► 3 scattering only in the examination of nuclear reactions is called democratic 
decay approximation. 

Nuclear structure studies within the framework of democratic decay approxima- 
tion have been started in 70th by Prof. R.I.Jibuti with collaborators. The reviews 
of the results obtained by the Tbilisi group can be found in refs. ^'^. Recently a 
considerable progress has been made in both theoretical and experimental studies 
of democratic decays. Below we shall briefly review some of the recent results. In 
particular, we shall discuss the results of experimental studies of democratic decays 
oi A = 6 nuclei via the three-body channel a -|- A^ -|- A^, theoretical investigations of 
monopole excitations of ^^C nucleus in the cluster model a + a + a and ^'^Li structure 
studies in the cluster model ^Li-|-n -|- n. 

2. Democratic decay of A = 6 nuclei 

Detailed experimental study of the three-body decay of T = 1, J'^ = 0^ and 2"^ 
states in ^He, ^Li and ^Be nuclei has been performed by the group from Kurchatov 
Institute (see ^'^ and references therein). The democratic decay concept has been for 
the first time experimentally tested in this investigation, and the democratic nature of 
the decay of these states via the only open channel a + N+N has been unambiguously 
demonstrated. 

The wave function in the democratic decay approximation has the following asymp- 
totics: 
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where the first term is a six-dimensional plane wave while the second one is an ingoing 
six-dimensional spherical wave; F is the 3-^3 scattering amplitude; Vab and rab,c are 
Jacobi coordinates, and g„^ and Qatc ^^^ associated canonically conjugated momenta; 
hyper-radius p = {rl^ + r\^^!'^ , and k = {ql^ + qlf^J^^^ = {2mEf/^/h- E is the 
energy measured from the three-body breakup threshold; S and Ms are the total spin 



of the system and its projection while Xsms is the spin part of the wave function. It 
is natural to expand the decay amplitude F in hyperspherical harmonics Ir = ^klm 
and to couple the three-body orbital momentum L and the spin S into the total 
angular momentum J . As a result, for the decay amplitude Fjm of a resonant state 
with the total spin J and its projection M we have the following expression ^'^: 

FjM = 22fKLS' ^KLMlXSMs\jj^j ■ (2) 

r 

In eq. (§), K is hypermomentum, the multi-index F = {K,L,Mi, S,'-f} plays the 
role of a channel index for the 3-^3 scattering, and the multi-index 7 stands for 
the set of all additional quantum numbers labelling different channels with the same 
quantum numbers K, L, Ml, S. It is convenient to use the Yamanouchi symbol r and 
the Young scheme [/] for 7 if there are three identical particles in the final state; 
otherwise orbital momenta lab and lab,c corresponding to Jacobi coordinates Vab and 
rab,c are usually used for 7. 

Eq. (^ is a generalization of the usual angular momentum decomposition of a 
two-body scattering amplitude and can be formally applied to any three-body decay. 
Note, that hypermomentum K is not an integral of motion, and due to the two-body 
interactions components labelled by different values of K mix together in the wave 
function, so, partial amplitudes /r = Jkls with various values of K are present in the 
expansion (^). For a two-body decay channel the expansion (§) is non-convergent; in 
the case when two-body subsystems of the decaying three-body system have sharp 
resonances, the convergence can appear to be very slow. If two-body subsystems of 
the three-body system have neither bound nor sharp resonant states, than the 3-^3 
scattering amplitude Fjm converges rapidly and only few terms play an important 
role in eq. (0). Thus only in the case of democratic decay of such system it is possible 
to extract partial amplitudes /r = Jkls fro^i the experimental data, i.e., to solve the 
problem that is somewhat similar to the problem of the spin assignment to a resonant 
state by the analysis of two-body scattering data. 

The hyperspherical harmonics Y^lml ^^^ functions of five angles in the formal six- 
dimensional momentum space {Qab^Qab^c}- 0^6 of the angles, ^d = a.Ycta.n{qab/qab,c), 
rules the energy distribution between two Jacobi coordinates. Therefore, in order 
to evaluate partial amplitudes Jkls^ ^^^ should analyze energy spectra of emitted 
particles. The analysis of ref. ^ of the three-body decay ^Be(0"^) -^ a + p + p will 
serve us as an example. 

^Be levels have been populated in the reaction ^Li(^He,^H)^Be at Eajf^ = 40 MeV. 
^H has been detected in the coincidence with a and/or p from the decay ^Be — *> 
a + p + p. The detection of ^H of appropriate energy made it possible to select 
the levels of interest. Taken from ref. ^ the spectrum of a-particles from the decay 
of the ^Be 0^ state in the beryllium rest frame is depicted on fig. 1. The peculiar 
form of the spectrum cannot be reproduced in the model of independent escape of 
the decay products (phase volume calculations), nor can it be reproduced in the 
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Figure 1: The a-particle spectrum (histogram) from ^Be(0"'") decay in the berylhum 
rest frame. The sohd hue is the hyperspherical harmonics description of the data. 
The figure is taken from ref. ^. 



model of sequential decay via the resonant states of p-p or a-p subsystems. At the 
same time, as is seen from fig. 1, the data are perfectly reproduced using expansion 
(H) and allowing only for the terms with K = and 2. Nevertheless, the set of 
the partial amplitudes {/^ls} cannot be uniquely determined from the analysis of 
the a-spectrum because various sets of {/kls} &^^ ^^^ ^^ ^^ ^^e data of the same 
quality. One should involve in the analysis the energy spectrum of protons from 
the same reaction. The predictions for the proton spectrum obtained using various 
sets of {fJcLs} fitted to the a-spectrum, differ significantly. From the fit to both a- 
and p-spectra, the relative weights of (6 ± 5)%; (44 ± 12)%; and (50 ± 17)% for the 
components K = 0; K = 2, 5* = 0; and K = 2, S = 1, respectively, have been 
obtained in ref. ^. 

The set of partial amplitudes {/^lsI ^^^ ^e used for the predictions of the kine- 
matically complete experiment based on the triple-coincidence ^H +p + a measure- 
ments. As is seen from fig. 2, the set of {/kls} fitted to a- and p-spectra, perfectly 
matches the triple-coincidence data. Note, that the kinematical conditions of the 
triple-coincidence experiment ^ have been chosen in such a way, that the contribu- 
tion of the S = 1 decay has been strongly suppressed. Therefore, the fine agreement 
between the theoretical curve and the experimental data on fig. 2 confirms the good 
quality of the determination of the relative weights of i^ = and K = 2, S = 
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Figure 2: The histogram is the triple coincidence distribution measured in the kine- 
matically complete experiment; the solid line is the prediction based on the a and 
proton spectra calculated using hyperspherical harmonics expansion. The figure is 
taken from ref. *^. 

components only. 



3. Democratic decay of 0^ states in ^^C 

In the previous section we have discussed the application of the democratic decay 
approximation to the analysis of experimental data. Now we shall discuss the use of 
the true three-body scattering theory in calculations. 

It is well known that the giant monopole resonance generates a well-pronounced 
peak of the monopole excitation strength in heavy nuclei exhausting an essential 
fraction of the energy- weighted sum rule (EWSR). In light nuclei, the giant monopole 
resonance is very weak and exhausts, e.g., in ^^C, only ~ 5% of EWSR '^. However, 
the strong peak of £"0 strength exhausting ~ 90% of EWSR has been predicted in 
various microscopic calculations (see, e.g., ref. ®). It has been shown in ref. ^ that 
this discrepancy should be attributed to the continuum spectrum effects, i.e., to the 
decay properties of the giant monopole resonance. 

The 3q; cluster model has been used in ref. ^ . The ground state wave function has 
been obtained by diagonalization of the three-body Hamiltonian using the harmonic 
oscillator basis. For 0"*" continuum spectrum states only the decay channel ^^C -^ 



a + a + a has been allowed for. This channel is the single open channel for low-lying 
0+ states and is supposed to dominate for 0"*" states with excitation energies of up 
to ~ 20 MeV. Continuum spectrum wave functions have been calculated by means 
of the oscillator representation of scattering theory ( J-matrix method) ^° which had 
been generahzed on the case of the true few-body scatting in ref. ^^. 

In the J-matrix approach, the continuum spectrum wave function ipl^ for a given 
three-body decay channel F = {K, 7}| is expressed as an infinite series expansion in 
six-dimensional harmonic oscillator functions (yC„x, 

oo 

^1 = E Dl^{E)v^K. (3) 

n=0 

The basic approximation of the J-matrix approach is in neglecting potential en- 
ergy matrix elements V^-J^ = (V12 -l- V13 -l- ^23)^^ with a large number of oscillator 
quanta N = 2n + K > N and/or A^' = 2n' + K' > N; the infinite kinetic energy 
matrix T is not truncated. The exact wave functions of the Hamiltonian with the 
infinite kinetic energy matrix and truncated potential energy matrix for any positive 
energy E can be found by algebraic methods. To do it, one should first diagonalize 
the truncated Hamiltonian matrix, i.e., calculate its eigenvalues Ex and corresponding 
eigenvectors C^p. Next the matrix elements V^ p (E) of the matrix V, 

K'f{E)=j:^^^, (4) 

^ Ex - E 

are calculated. The 3^3 scattering ^-matrix can be easily computed by the expres- 
sion 

S = {A^+Y^A^-^ , (5) 

where matrix elements of matrices A^~^^ and A^~^ , 

4±) = 5pp, ci%{E) + V:'f{E) Tlly^, Cltl^K'iE) ; (6) 

u = |(iV — K) and z/' = |(iV — K') are the truncation boundaries of the Hamiltonian 
matrix in the channels F = {K,'~f} and F' = {i^',7'}, respectively; the kinetic energy 
matrix elements T^^^^^ = —\hoJ J {n + l){n + K + 3); the eigenvectors C\^j^{E) and 

C!^j({E) of the infinite kinetic energy matrix T can be easily calculated by analytical 
expressions or by recurrent relations (see ref. ^^ for details). 

Finally, for the wave function with outgoing six- dimensional spherical wave in 
the channel F' and ingoing six- dimensional spherical waves in all allowed channels F, 

^Hereafter we shall omit the spin part of the wave functions and quantum numbers S, L, and 
M]^ in the channel index F = {K,L,AI]^,S,'^} and in the indexes of the functions used. In the 
case of spinless particles (e.g., a-particles), S — 0, L, and M^ are integrals of motion, therefore, all 
functions are labelled with the same values of S, L, and M^ and channels T = {K, L, M^, S, 7} and 
r' = {K\ L', M'^,,S', 7'} with L^ L' and/or M ^ M'^, are not coupled. 
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the expansion coefficients Dl^^. (E) = Dlj^{E) [see eq. (|])] in the asymptotic region 
n> u are calculated by the expression 

D^lPiE) = ^[ 5rr' ci+i(i?) - Srr dJ{E) ] , (7) 

while for calculation of D^^ \E) in the interaction region n < v the expression 

D^nv\E) = - E K^v"{E) Tj;;;,,,^, Df:i^,„{E) (S) 

r" 

is used. 

Ali-Bodmer aa potential ^^ has been used in ref. ^; Coulomb interaction has been 
neglected. The truncation boundary iV = 18 has been used in the calculations, i.e., 
the truncated Hamiltonian matrix accounted for all 3a oscillator states up to I'^huo. 
The convergence for both the ground and continuum O"*" states is achieved at A^ ~ 18 
provided that huo = 10 MeV. In the asymptotic region n > z/ all channels F with 
K < 10 have been allowed for. However, in the interaction region n < z/ we have 
accounted for all oscillator states up to 18^ti;, thus the hyperspherical harmonics with 
i^ < 18 have been allowed for in the interaction region. The inclusion of components 
with 7^ > 10 in the channel space have not produced changes in the results. 

The results of the calculations of ref. ^ are presented on fig. 3. The solid line in 
fig. 3 represents the ii^-dependence of the £^0 transition strength F[E) = E \Mq\ joao 
to the states above the ^^C — > 3a threshold. Here, Mq is the matrix element of the 
EO transition from the ground state to the states with energy E. The dash-dotted 
curve displays the sum of the contributions to EWSR from all states with excitation 
energies less than E^ i.e., the function 

S{E) = ^t^ X 100% , (9) 

0"oo 

where 

a{E)= [ E'\MofdE' , (10) 

Jo 

the monopole sum rule Uoo = cr(oo) = ^^ (-R^)o; and ((-R^)q) is the rms radius of 
the ^^C nucleus. The dashed step-like line is the plot of the function S{E) calculated 
disregarding the effect of the 3a continuum. In this case, the height of each step is 
equal to the contribution to the EWSR from the corresponding quasistationary state 
obtained by the diagonalization of the truncated Hamiltonian matrix. It is seen that 
the most essential contributions to EWSR are connected with the quasistationary 
levels with excitation energies 6.11, 10.39, 13.60 and 22.30 MeV (31%, 29%, 13% and 
13%, respectively). All the remaining levels with excitation energies up to 35 MeV 
exhaust about 7% of the EWSR. 

The ^^C binding energy in calculations of ref. ^ has been found to be Eb = 
7.27 MeV (the experimental value is 7.27 MeV). The ffist excited O'*' state with the 
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Figure 3: The solid line displays the monopole strength function F{E) = E \Mq\ JGoq 
for the states in continuum. The dash-dotted line is a plot of the function S{E) [see 
eq. @], i.e., it displays the fraction of monopole EWSR exhausted by the states 
with excitation energies less than E. The dashed line is the same, but obtained by 
neglecting the effect of continuum. The experimental energies and the widths of the 
£^0 resonances are indicated by crosses ^ and circles ^^. The arrow indicates the ^^C 
-^ So; reaction threshold. The figure is taken from ref. ^. 



excitation energy Ex = 6.11 MeV belongs to the discrete spectrum, so, the continuum 
does not affect its properties. The experimental value for Ei is 7.66 MeV and this 
state belongs to the 3a continuum. However, its width is very small (F = 9 eV) and 
its properties are rather close to those of the discrete spectrum states. 

The allowance for the continuum slightly changes the properties of the 10.39 MeV 
0^ state. Namely, the F{E) maximum shifts down to £^ = 9.9 MeV, the state gets 
some width and its contribution to EWSR decreases to 26%. 

However, in the 13-25 MeV excitation region the effect of the continuum appears 
to be essential. The levels in this excitation energy region get considerable widths and 
merge each other, forming very broad (of ~ 12 MeV in width) plateau of a peculiar 
shape, i.e., the strength function F{E) is almost constant between 13 and 25 MeV 
except for a bump of the width of 2-3 MeV at an excitation energy of about 20 MeV. 
The bump should be related to the 20 MeV monopole resonance observed in refs. ^^'^. 
33% of the EWSR is exhausted in the 12-26 MeV excitation energy region, with 16% 



of the EWSR being exhausted in the 18-22 MeV interval. 

The approach, as is seen from fig. 3, reproduces quite correctly the positions of 
the strong EO absorption regions. The broad and smooth pedestal of the 20 MeV EO 
excitation strength maximum can be easily misinterpreted in experimental studies as 
a background. The shape of the background that has been subtracted by Eyrich et 
al ^ in their study of the 20 MeV resonance, is similar to the shape of the pedestal on 
fig. 3. If the experimental background is a real manifestation of the smooth pedestal 
of the EO strength, than it should be higher for small scattering angles. From fig. 3 
of ref. ^ it is seen, that the background is really enhanced at small angles. So, the 
experimental spectra of ref. ^ is in qualitative agreement with theoretical results of 
ref.^. The further investigation of the nature of the background is needed. 

The results presented on fig. 3 give the explanation why the giant monopole res- 
onance is nearly not seen in ^^C There are two physical reasons. First one is the 
strong fragmentation of the EO strength in ^^C As a result, more than 50% of the 
EWSR is exhausted by the 7.65 MeV and 10 Mev resonances, and not too much of 
the EO strength is left for the giant monopole resonance. The second reason is the 
continuum spectrum effect. Due to the democratic decay, the giant monopole reso- 
nance gets a large width and merge with other O"*" states. As a result, the £^0 strength 
becomes very smooth, and the giant monopole resonance manifests itself only as a 
non-prominent bump on the smooth pedestal. About 15% of the EWSR is exhausted 
in the vicinity of 20 MeV bump, but if the pedestal is treated as a background, than 
the contribution of the resonance to the EWSR is severely underestimated. 

4. ^^Li structure and the soft dipole mode 

Light neutron-rich nuclei like ^He, ^He, ^^Li, ^^Be, etc., are extensively studied 
now (see recent review papers ^^'^^ and references therein). These nuclei are often 
called exotic nuclei; they are also often called Barromean nuclei, because they have 
a well-pronounced cluster structure core + n + n and after removing a neutron from 
any of these nuclei, the remaining system core + n becomes unbound.^ Two excess 
loosely-bound neutrons have a wide spatial distribution with long tail reaching far out 
and form the so-called neutron halo. The neutron halo manifests itself in a number 
of unusual properties of exotic nuclei, in particular, in a large value of the rms radius 
of the nucleus, in the existence of the soft dipole mode, etc. 

It is natural to make use of the democratic decay approximation in the study of 
Barromean nuclei. In ref. ^^, ^^Li has been studied in the cluster model ^Li + n + n. 
It has been shown in ref. ^^ that the democratic decay approximation within the 
framework of the J- matrix approach described in the previous section, is advantageous 
for the investigations not only of the continuum states, but of the ground state as 



■I- The heraldic symbol of the dukes of Barromeo is three rings connected in such a manner that 
after removing of any of the rings the two others become disconnected. 



well. 

In a usual variational calculation with the oscillator basis, it is impossible to 
describe the slowly-dying asymptotic tail of the halo neutron wave functions. In 
order to describe the neutron halo properties of the ^^Li nucleus, one is pushed to 
account for the oscillator states with a very large number of oscillator quanta. In the 
J-matrix approach, the infinite number of the oscillator states can be allowed for. 
For a bound state, the expansion coefficients Dnr{E) = D^j^{E) [see eq. (|)] in the 
asymptotic region n>v are of the form^^: 

Dr^E) = ar cl^kE) . (11) 

The bound state asymptotic normalization constants, ar, are calculated by handling 
numerically the set of linear homogeneous equations^^ 

E 4r' «r' = . (12) 



detAW = 0. (13) 



Eq. (|T2D can be solved if only 



Using eq. (0) the energy of a bound state can be easily calculated. This method of 
calculation of the bound state energies is equivalent to the location of the S'-matrix 
pole^^ corresponding to the bound state, as is easily seen from eq. (|). 

The expansion coefficients Dnr{E) = D^p^[E) for the inner region n < v are 
obtained using eq. (P). Finally, the wave function should be normalized. 

The calculations^^ have been performed using the shallow Gaussian n-^Li and 
n-n potentials of ref. ^^. Only the states with the spin of the halo-neutron pair 5* = 
have been considered. The three-body angular momentum of the ground state of ^^Li 
have been supposed to be equal to L = 0. The parameter hu = 7.1 MeV has been 
used in the calculations. This value corresponds approximately to the minimum of 
the ground state energy Eq. 

The results for the ^^Li ground state for various values of the truncation parameter 
N are presented in the table 1. The variational ground state energies, Eq , obtained 
by the pure diagonalization of the truncated Hamiltonian matrix are listed in the 
second column, while the J-matrix results, Eq, which are the solutions of the eq. (pTSD 
are listed in the third column. The values of the ^^Li r.m.s. radius, < r^ >^( and 
< r^ >i( , obtained by the pure diagonalization of the truncated Hamiltonian matrix 
and with the allowance for the asymptotic tail of the wave function, respectively, are 
presented in the 4-th and the 5-th columns of the table. It is seen, that by locating the 
S'-matrix pole using eq. ([T3|) that is equivalent to the allowance for the long asymptotic 
tail of the wave function, the convergence of the ground state is essentially improved. 

The results presented in the table 1 have been obtained using Lanczos smoothing 
of the three-body potential energy matrix^^. The smoothing improves the convergency 
of calculations of the scattering characteristics. Nevertheless, the smoothing causes 
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Table 1: ^^Li ground state properties (see text for details). 





Ground state energy, 


Neutron halo 


Truncation 


MeV 


mean square radius 


boundary A^ 




< r^ >li\ fm 


^0 


Eo 


< r^ >[{' ^'^ 


< ^ >ii 


12 


-0.012 


-0.150 


2.83 


3.31 


16 


-0.116 


-0.199 


2.91 


3.29 


20 


-0.171 


-0.225 


2.98 


3.31 


24 


-0.202 


-0.240 


3.04 


3.32 


Experiment 


-0.247±0.080 


3.16±0.11 



Table 2: The effect of the Lanczos smoothing on the ground state energy of ^^Li (see 
text for details). 



N 


10 


12 


14 


16 


18 


20 


22 


24 


Eo 


-0.119 


-0.150 


-0.181 


-0.199 


-0.215 


-0.225 




-0.240 


Eq 


-0.136 


-0.266 


-0.204 


-0.262 


-0.232 


-0.262 


-0.249 


-0.263 



the underestimation of the binding energy of a three-body system. The ^^Li ground 
state energy obtained in the J-matrix calculations without the use of the smoothing, 



E^ 



(w.s.) 



shows a staggering as A^ increases 



(see table 2). Eq'^ is seen to converge to 



the value of —0.263 MeV. The smoothing results in a usual-type monotonic decrease 
of Eq as N increases, but Eq''^ < Eq for all values of N. At the same time, the 
ground state wave functions obtained with the smoothing and without it, differ very 
slightly. So, below we shall discuss only the results of calculations with the smoothed 
potential energy matrix. 

The quality of the convergence of the wave function in the J-matrix approach 
can be demonstrated by the plots of momentum distributions. Figure 4 presents the 
transverse momentum distribution of the cluster ^Li in ^^Li that is currently supposed 
(the so-called Serber reaction mechanism) to be proportional to the ^Li transverse 
momentum distribution -^ in the fragmentation of high-energy ^^Li beams on target 
nuclei. The experimental data for ^^Li fragmentation^^ are also presented on the 
figure. The results of J-matrix calculations with A^ = 16 and A^ = 24 are so close to 
each other that their plots on fig. 4 coalesce and cannot be distinguished. 

It is seen that in calculation of the ground state, the allowance for the wave 
function asymptotics is very important for a weekly-bound system like ^^Li. The 
terms of expansion (0) with the number of oscillator quanta A^ ~ 100 that cannot be 
obtained in the usual oscillator-basis variational calculations, play an essential role 
in the formation of the transverse momentum distribution, r.m.s. radius, etc. The 



§Note, that J-matrix calculation is not a variational one, so, the ground state energy may behave 
non-monotonically as N increases. 
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Figure 4: The ^Li transverse momentum distribution in the ground state of ^^Li. 1 
— variational calculations, 2 — J-matrix calculations. Experimental data for the 
^Li transverse momentum distribution in the fragmentation of ^^Li 800 MeV/nucleon 
beam are taken from ref. ^^. 

convergence of < r^ >^( , transverse momentum distribution and other properties 
of the wave function (e.g., of the weights of its components) in the full J-matrix 
calculation is very good. Nevertheless, it is seen that the r.m.s. radius converges 
to a value that is somewhat larger than the experimental one, and the calculated 
transverse momentum distribution is narrower than the experimental one. These 
shortcomings can be overcome by the adjustment of n-^Li potential. We have not 
aimed to fit the potential to the ^^Li properties, we have just take its parameters 
from ref. ^^. The differences between the calculated and experimental values of r.m.s. 
radius arise from the fact that in ref. ^^ the potential has been fitted to the ^^Li 
r.m.s. radius and the binding energy in the variational calculation that is unable to 
reproduce these quantities with high accuracy. 

The reduced probability of the El transition obtained in the cluster model, 
B{E1; Ef — Eq), is displayed on the figure 5. In the J-matrix approach, zero-width 
pure variational peaks of B{E1; Ef — Eq) gain finite widths and shift to lower en- 
ergies. It is seen that in the calculation of B{E1; Ef — Eq) it is important to allow 
for the asymptotics not only for continuum states, but for the ground state as well. 
Due to the wide spatial distribution of the halo neutrons density, large distances con- 
tribute essentially to the El strength in the low-energy region. To account for this 
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Figure 5: Cluster B{E1; g.s. -^ continuum) in ^^Li. Zero-width peaks with diamonds 
on the top (1) given in arbitrary units, have been obtained in the pure variational 
approach for both the ground and exited states, i.e. the wave functions for the 
initial and the final states have been obtained by the diagonalization of the truncated 
Hanailtonian matrix. The curve with dots (2) has been calculated with the variational 
ground state wave function while the continuum spectrum wave function has been 
obtained by the J-matrix approach. The solid curve (3) is the result of J-matrix 
calculations for both the ground and continuum spectrum states. 

contribution and to get a convergence, one should allow for the ground and contin- 
uum state wave function components with very large values of the total number of 
oscillator quanta A^ ~ 2000^] in calculations of B{E1; Ef — Eq) . The large-distance 
El transitions enhance the peak of the cluster B{E1; Ef — Eq) and shift it to a lower 
energy. Obviously, this peak should be associated with the soft dipole mode. Thus, 
the neutron halo manifests itself in the appearance of the soft dipole mode that arises 
from the shift and enhancement of the B{E1] Ef — Eq) peak. This effect appears 
to be very important in the calculation of the ^^Li electromagnetic dissociation cross 
section (see below). 

Large distances contribute essentially to E2 and EO transitions, too. The shift 
and the enhancement of B{EO] Ef — Eq) and of B{E2; Ef — Eq) are even more 
pronounced. Nevertheless, E2 and EO transitions do not play an important role in 
the electromagnetic dissociation, and we shall not discuss them here. 

Figure 6 shows the comparison of the results of our calculations^^ of the clus- 

^The classical turning point for the basis function with the total number of oscillator quanta 
N = 1000 is at a distance of «108 fm from the origin. 
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Figure 6: Comparison of our results^^ for B{El;g.s. -^ continuum) in ^^Li with 
results of other authors. 1 — J-matrix method^^; 2 — ref. ^^; 3 — ref. ^^; 4 — 
experimental data parametrization of ref. ^°. 

ter B{E1; Ef — Eq) with the parametrization of experimental data of ref. ^°. The 
agreement is reasonable. The form of the B{E1; Ef — Eq) peak is well reproduced, 
the discrepancy in the position of the B{E1; Ef — Eq) maximum is supposed to be 
eliminated by the adjustment of the potentials. The results of the B{E1; Ef — Eq) 
calculations of refs. ^^'^^ are also depicted on fig. 6. The three-body cluster calcu- 
lations with the allowance for democratic decay channels of ref. ^^ performed using 
usual coordinate-space hyperspherical harmonics approach, nicely reproduce the en- 
ergy of the soft dipole mode. Nevertheless, the form of the peak in our calculations is 
reproduced better. Note, that the authors of ref. ^^ used another set of the potential 
parameters. The calculations of ref. ^^ with the allowance for two-body decay chan- 
nels only, failed to reproduce the form of the B{E1; Ef — EQ) peak and underestimate 
the energy of the soft dipole mode.. 

The soft dipole mode exhausts about 90% of the cluster EWSR, a^""**, associated 
with the excitation of the cluster degrees of freedom only. Nevertheless, it is easy to 
show that 

1.3%, (14) 



oo 

^"12 



where o"^* is the total EWSR accounting for excitations of all nucleons. So, the 
contribution from the soft dipole mode to the total EWSR is relatively small. In the 
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vicinity of the B{E1; Ef — Eq) maximum at the excitation energy E k, 1 MeV only 
~8% of the total EWSR is exhausted. Nevertheless, the account for the soft dipole 
mode results in an essential increase of the electromagnetic dissociation cross section 
of 0.8 GeV/nucleon ^^Li beams on Pb and Cu targets. The wide space distribution of 
the halo neutrons density is also well-manifested in the electromagnetic dissociation of 
^^Li beam. The contribution of the large-distance i?A-transitions to the cross section is 
about 50%. The calculations have been performed by the equivalent photon method. 
The only parameter of the method is the minimal value of the impact parameter hmin- 
We use for hmin the values of 9.0 fm for Pb and 6.8 fm for Cu target nuclei, respectively. 
These quantities are the sums of the ^^Li and target nucleus charge radii. With these 
values of hmin we obtain for the electromagnetic dissociation cross sections the values 
of 0.966 barn for the Pb target and 0.132 barn for the Cu target; the corresponding 
experimental values are 0.890 ±0.110 barn and 0.21±0.04 barn, respectively^^. EO 
and E2 transitions give only 1.2% contribution in the cross sections. 
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